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We study a novel configuration for displacement detection consisting of a nanomechanical res- 
onator coupled to both, a radio frequency superconducting interference device (RF SQUID) and to 
a superconducting stripline resonator. We employ an adiabatic approximation and rotating wave 
approximation and calculate the displacement sensitivity. We study the performance of such a 
displacement detector when the stripline resonator is driven into a region of nonlinear oscillations. 
In this region the system exhibits noise squeezing in the output signal when homodyne detection 
is employed for readout. We show that displacement sensitivity of the device in this region may 
exceed the upper bound imposed upon the sensitivity when operating in the linear region. On the 
other hand, we find that the high displacement sensitivity is accompanied by a slowing down of the 
response of the system, resulting in a limited bandwidth. 
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I. INTRODUCTION 

Resonant detection is a widely employed technique in a 
variety of applications. A detector belonging to this class 
typically consists of a resonator, which is characterized 
by a resonance frequency luq and characteristic damping 
rates. Detection is achieved by coupling the measured 
physical parameter of interest, denoted as a;, to the res- 
onator in such a way that becomes effectively x depen- 
dent, that is luq — loq (x). In such a configuration x can be 
measured by externally driving the resonator, and mon- 
itoring its response as a function of time by measuring 
some output signal X (t) . Such a scheme allows a sensi- 
tive measurement of the parameter x, provided that the 
average value of X (t), which is denoted as Xq, strongly 
depends on luq, and provided that wq, in turn, strongly 
depends on x. These dependencies are characterized by 
the responsivity factors R = \dXo/dujo\ and \dujo/dx\ re- 
spectively. Resonant detection has been employed before 
for mass detection jlj , quantum state readout of a quan- 
tum bit [1, d, 0, 13| , detection of gravitational waves , 
and in many other applications. 

In general, any detection scheme employed for moni- 
toring the parameter of interest x can be characterized by 
two important figures of merit. The first is the minimum 
detectable change in x, denoted as 6x. This parameter 
is determined by the above mentioned responsivity fac- 
tors, the noise level, which is usually characterized by 
the spectral density of X (t) , and by the averaging time 
T employed for measuring the output signal X (t) . The 
second figure of merit is the ring-down time ^rdj which 
is a measure of the detector's response time to a sudden 
change in x. 

In general, the minimum detectable change dx is pro- 
portional to the square root of the available bandwidth, 

1/2 

namely, to {2it/t) ' . It is thus convenient to charac- 



terize the sensitivity by the minimum detectable change 
Sx per square root of bandwidth, which is given by 

1/2 

= Sx/ {2Tr /tY'. Under some conditions, which will 
be discussed below in detail, the smallest possible value 
of Pa; is given by P, 0] 
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where kgT is the thermal energy, Uq is the energy stored 
in the resonator, and Q is the quality factor of the res- 
onator. One of the assumptions, which are made in order 
to derive Eq. ([1]), is that the response of the resonator 
is linear. We therefore refer to the value of Px given by 
Eq. Il]) as the standard linear limit (SLL) of resonant 
detection. Under the same conditions and assumptions, 
the ring-down time is given by 



^RD = — 

OJo 



(2) 



As can be seen from Eq. ([T]), sensitivity enhancement 
can be achieved by increasing Q, however, this unavoid- 
ably will be accompanied by an undesirable increase in 
the ring-down time (see Eq. ([2])), namely, slowing down 
the response of the system to changes in x. Moreover, Eq. 
(H]) apparently suggests that unlimited reduction in P^ 
can be achieved by increasing Uq by means of increasing 
the drive amplitude. Note however that Eq. ([1]), which 
was derived by assuming the case of linear response, is 
not applicable in the nonlinear region. Thus, in order 
to characterize the performance of the system when non- 
linear oscillations are excited by an intense drive, one 
has to generalize the a naly sis by taking nonlinearity into 
account [3, iZ, S i, [M UlI • 

In the present paper we theoretically study a novel 
configuration for resonant detection of displacement of a 
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nanomechanical resonator [l^, [l3|, which is coupled to 
both, a radio frequency superconducting interference de- 
vice (RF SQUID) and to a superconducting striphne res- 
onator. Similar configurations have been studied recently 
in [ll, [H, [H, [13, [il, Ii9]. Wc employ an adiabatic ap- 
proximation and a rotating wave approximation (RWA) 
to simplify the analysis and calculate the displacement 
sensitivity. We first consider the case where the response 
of the stripline resonator is linear and reproduce Eqs. ([1]) 
and ^ in this limit. However, we find that the response 
becomes nonlinear at a relatively low input drive power. 
Next, we show that Px in the nonlinear region may be- 
come significantly smaller than the value given by Eq. 
dU), exceeding thus the SLL imposed upon the sensitiv- 
ity when operating in the linear region. On the other 
hand, we find that the enhanced displacement sensitiv- 
ity is accompanied by a slowing down of the response of 
the system, resulting in a limited bandwidth, namely, a 
ring-down time much longer than the value given by Eq. 



II. THE DEVICE 

The device, which is schematically shown in Fig. 
[U consists of an RF SQUID inductively coupled to a 
stripline resonator. The loop of the RF SQUID, which 
has a self inductance A, is interrupted by a Josephson 
junction (JJ) having a critical current Ic and a capaci- 
tance Cj. A perpendicularly applied magnetic field pro- 
duces a flux $e threading the loop of the RF SQUID. The 
stripline resonator is made of two identical stripline sec- 
tions of length 1/2 each having inductance Lt and capac- 
itance Ct per unit length and characteristic impedance 
Zt = \/ Lt /Ct- The stripline sections are connected by 
a doubly clamped beam made of a narrow strip, which 
is freely suspended and allowed to oscillate mechanically. 
We assume the case where the fundamental mechanical 
mode vibrates in the plane of the figure and denote the 
amplitude of this flexural mode as x. Let m be the ef- 
fective mass of the fundamental mechanical mode, and 
ujm its angular resonance frequency. The suspended me- 
chanical beam is assumed to have an inductance Lb (in- 
dependent of x) and a negligible capacitance to ground. 
The coupling between the mechanical resonator and the 
rest of the system originates from the dependence of the 
mutual inductance M between the inductor and the 
RF SQUID loop on the mechanical displacement x, that 
is M = M {x). 



A. Transmission Line Resonator 

The transmission line is assumed to extend from y = 
—1/2 to y = 1/2, and the lumped inductance Lj, is located 
at ?; = 0. Consider the case where only the fundamental 
mode of the resonator is driven. Disregarding all other 
modes we express the voltage V (y, t) and current X (y, t) 




FIG. 1: (Color online) The device consists of a nanomechan- 
ical resonator coupled to both, a radio frequency supercon- 
ducting interference device (RF SQUID) and to a supercon- 
ducting stripline resonator. 



along the transmission line as 
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(4) 



where Lp represents the flux in the lumped inductor at y = 
0. The value of k, is determined by applying Faraday's 
law to the lumped inductor at y = [18| 



td kI Lb 



(5) 



For the fundamental mode the solution is in the range 
TT < < 27r (see Fig. [2]). 



B. Inductive Coupling 

The total magnetic flux <f> threading the loop of the 
RF SQUID is given by 



(6) 



where $j represents the flux generated by both, the cir- 
culating current in the RF SQUID 1^ and by the current 
in the suspended mechanical beam If, 



^i=IsA + MIb 



(7) 



where A is the self inductance of the loop. Similarly, the 
magnetic flux ip in the inductor Lb is given by 



(f = IbLb + Mis ■ 
Inverting these relations yields 

Lb'i'i-Mip 



h = 



kLb{l-alj) 



(8) 

(9) 
(10) 
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where 



M 



(11) 



The gauge invariant phase across the Josephson junc- 
tion 6 is given by 

27r$ 

= 2nn~—, (12) 
q?o 

where n is an integer and $0 = h/2e is the flux quantum. 
We set n = 0, since observable quantities do not depend 
on n. 



C. Capacitive and Inductive Energies 

Assuming that the only excited mode is the funda- 
mental one, the capacitive energy stored in the stripline 
resonator is found using Eqs. ([3]) and (O 



^ J-l/2 

where Ce, which is given by 



(13) 



(14) 



represents the effective capacitance of the stripline res- 
onator. The factor i?, which is defined by 



sin {kI) 
kI + sin (kI) 



(15) 



can be calculated by solving numerically Eq. (O (see 
Fig. [2]). Similarly, the inductive energy stored in the 
resonator, excluding the energy stored in the lumped in- 
ductor Lb at y ^ 0, is found using Eqs. ([4]) and ^ 



2 J-i/2 2Le 



where Le, which is given by 
1 

1, 



(16) 



(17) 



represents the effective inductance of the stripline res- 
onator excluding the lumped element at y = 0, and where 

(18) 



V LtCt 



The inductive energy stored in the RF SQUID loop 
and the lumped inductor Lb is calculated using Eqs. © 
and HOl) 



2Lb 



A M 

M Lb 



2A(l-0 



(19) 
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FIG. 2: (color online) The factor nl (blue solid line) and ■& 
(green dashed line) calculated using Eq. ([5]) as a function 
Lt/Lrl. 



III. LAGRANGIAN AND HAMILTONIAN OF 
THE CLOSED SYSTEM 

Here we derive a Lagrangian for the closed system 
consisting of the nanomechanical resonator, stripline res- 
onator and the RF SQUID. The effect of damping will 
be later taken into account by introducing coupling to 
thermal baths. The Lagrangian of the closed system is 
expressed as a function of x, ip and $ and their time 
derivatives (denoted by overdot) 



C = 



2 2 2 

where the potential terms are given by 



-Uo-Ui, 



Uo 



? 2 
2 



^ <i>0/cCOS^ 



2A(1 



2tt 



(20) 



(21a) 



(21b) 



Using Eqs. (O, (O, pO| and (fTTjl . the corresponding 
Euler - Lagrange equations can be expressed as 



nit + mu}„-^x — I sib 
Lf. 

Cj$-t-/s -l-/csin- 



dAf 

da; 

27r$ 
$0 



, 
, 

. 



(22a) 
(22b) 

(22c) 



The interpretation of these equations of motion is 
straightforward. Eq. (|22a[) is Newton's 2nd law for the 
mechanical resonator, where the force is composed of the 
restoring elastic force —muj^x and the term due to the 
dependence of M on x. Eq. (I22bp relates the current 
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in the suspended beam with the currents in the ef- 
fective capacitor Ce and the effective inductor L^. Eq. 
(|22c|) states that the circulating current Ig equals the 
sum of the current Ic sin 9 through the J J and the cur- 
rent CjVj through the capacitor, where the voltage Vj 
across the JJ is given by the second Josephson equation 
l/j = ($o/27r) ^. 

The variables canonically conjugate to x, Lp and $ are 
given hy p = mi, q = Ce^ and Q — Cj<f> respectively. 
The Hamiltonian is given by 



The total wave function is expanded as 



where 



9 ? 



Hi = 



2Cj 



Ui 



(23) 



(24) 



(25) 



Quantization is achieved by regarding the variables 
{x,p,(p,q,^,Q} as Hermitian operators satisfying Bose 
commutation relations. 



IV. ADIABATIC APPROXIMATION 

As a basis for expanding the general solution we use 
the eigenvectors of the following Schrodinger equation 



Hi \n {x, ip)) = £„ {x, ip) \n {x, p)) 



(26) 



where x and Lp are treated here as parameters (rather 
than degrees of freedom). The local eigen- vectors are 
assumed to be orthonormal 



{m{x,ip) \n{x,ip)) = (S„ 



(27) 



The eigenenergies e„ (x, (p) and the associated wave- 
functions ipn are found by solving the following 
Schrodinger equation 



-/3c 



92 



where 



^0 



2/3l cos( 



Eo 



27r$ 
$0 

27r$e 

$0 
27rA/c 

$0 

CjE() 
87r2A 



2ttM (x) ip 



(28) 

(29) 

(30) 

(31) 
(32) 

(33) 
(34) 



V' = X!^" (x,(^,t) \n) 



(35) 



In the adiabatic approximation [20| the time evolution 
of the coefHcients ^„ is governed by the following set of 
decoupled equations of motion 



[Ho + En {X, (f)] ^„ = ihi„ 



(36) 



Note that in the present case the geometrical vector po- 
tential [20[ vanishes since the wavefunctions ipn {4>) can 
be chosen to be real. The validity of the adiabatic ap- 
proximation will be discussed below. 



V. TWO LEVEL APPROXIMATION 

In what follows we focus on the case where \(f>o\ ^ 1 
and Pl (1 — ck^) > 1. In this case the adiabatic poten- 
tial u ((/)) given by Eq. ([^ contains two wells separated 
by a barrier near = 0. At low temperatures only the 
two lowest energy levels contribute. In this limit the lo- 
cal Hamiltonian Hi can be expressed in the basis of the 
states \~r\) and \rx), representing localized states in the 
left and right well respectively having opposite circulat- 
ing currents. In this basis, Hi can be expressed using 
Pauli's matrices 



(37) 



The real parameters 77 and A can be determined by solv- 
ing numerically the Schrodinger equation (|28p 14]. The 
eigenvectors and eigenenergies are denoted as Hi\±) = 
e± |±), where 



(38) 



VI. ROTATING WAVE APPROXIMATION 



Consider the case where $e — $o/2, that is 



yo = 



2'kM [x) ip 



(39) 



and assume that adiabaticity holds and that the RF 
SQUID remains in its lowest energy state. In this case, as 
can be seen from Eq. (|38p . expansion of e_ (x, ip) yields 
only even powers of tp. These even powers of ip can be 
expressed in terms of the annihilation operator 



UJeip - 



\/Ce 



(40) 



and its Hermitian conjugate At, yielding terms oscillat- 
ing at frequencies 2nWe, where n is integer. In the RWA 
such terms are neglected unless n = since the effect 
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of the oscillating terms on the dynamics on a time scale 
much lon ger than a typical oscillation period is negligi- 
bly small [2l[. Moreover, constant terms in the Hamilto- 
nian are disregarded since they only give rise to a global 
phase factor. Displacement detection is performed by 
externally driving the fundamental mode of the stripline 
resonator. To study the effect of nonlinearity to lowest 
order we keep terms up to 4th order in ip. On the other 
hand, since the mechanical displacement is assumed to 
be very small, we keep terms up to 1st order only in x. 
Thus, in the RWA the Hamiltonian Tip + £- is given by 

T^RWA = huJr„.Nrn + ^Wq {x) + HK NI , (41) 

where and N^, = A\Af, are number operators of the 
mechanical and stripline resonators respectively. 



K ~ — — 



^ / „dlogM 

LUO (X) = CJe — >'h [ i- + 2 ; X 



dx 



^2 = ^ 



A T] 2ttMo / h 



h \A $oLb V 2CeCj, 



(42) 
(43) 
(44) 



A. Equation of Motion 

Using the standard method of Gardiner and CoUett 
[23| . and applying a transformation to a reference frame 
rotating at angular frequency ujp 

= Ce-^'^p* , (46) 

yield the following equation for the operator C 

dC 



dt 



+ e = Fit) 



(47) 



where 



and Mo = A/(0). 



e (c, ct) ^[j + i{l,o- ujp) + {iK + 73) etc] c 

(48) 

7 — 71+72, 4>i is the (real) phase shift of transmis- 
sion from the feedline into the resonator, and F it) is a 
noise term, having a vanishing average {F (t)) = 0, and 
an autocorrelation function, which is determined by as- 
suming that the three semi-infinite transmission lines are 
at thermal equilibrium at temperatures Ti, T2 and Ta 
respectively. 



VII. HOMODYNE DETECTION 



B. Linearization 



The stripline resonator is weakly coupled (with a cou- 
pling constant 71) to a semi infinite feedline, which guides 
the input and output RF signals. To model the effect of 
dissipation (both linear and nonlinear), we add two ficti- 
tious semi infinite transmission lines to the model, which 
allow energy escape from the resonator. The first trans- 
mission line is linearly coupled to the resonator with a 
coupling constant 72, and the second one is nonlinear ly 
coupled with a coupling constant 73 [23 |. 

The dependence of on x can be exploited for dis- 
placement detection. This is achieved by exciting the 
fundamental mode of the stripline resonator by launch- 
ing into the feedline a monochromatic input pump signal 
having a real amplitude bin and an angular frequency LUp 
close to the resonance frequency ujq. The output signal 
Cout reflected off the resonator is measured using homo- 
dyne detection, which is performed by employing a bal- 
ance mixing with an intense local oscillator having the 
same frequency as the pump frequency LOp, and an ad- 
justable phase That is, the normalized (with re- 
spect to the amplitude of the local oscillator) output sig- 
nal of the homodyne detector is given by 



Xd,, n. — 



(45) 



To proceed, we employ below some results of Ref. [22l |. 
which has studied a similar case of homodyne detection 
of a driven nonlinear resonator. 



Let C — B + c, where _B is a complex number for which 

e{B,B*)^0, (49) 

namely, B is a steady state solution of Eq. ([Tf|) for the 
noiseless case F — 0. When the noise term F can be 
considered as small, one can find an equation of motion 
for the fluctuation around B by linearizing Eq. (I47p 



— + Wc+Vc'=F, 
dt 



(50) 



where 



W = 



and 



59 

dC 



= -/ + i{uJo~LUp) + 2 {iK + 73) B*B 



C=B 



V = 



de 



9Ct 



= (iK + 73) B^ 



(51) 
(52) 



C=B 



C. Onset of Bistability Point 

In general, for any fixed value of the driving amplitude 
bin, Eq. can be expressed as a relation between \B\^ 
and LUp. When bi„ is sufficiently large the response of the 
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resonator becomes bistable, that is |B| becomes a multi- 
valued function of ujp. The onset of bistability point is 
defined as the point for which 



d\B 



^- = 

2 " 



= . 



d \B\ 



(53) 
(54) 



Such a point occurs only if the nonlinear damping is suf- 
ficiently small [1^ , namely, only when the following con- 
dition holds 



\K\ > VS' 



73 



(55) 



At the onset of bistability point the drive frequency and 
amplitude are given by 



473|i^| + 73(g^+7|) 
~ 37| 



(56) 



where W is the real part of W. Thus one has 



Xo.i^W 1± Jl 



m 



lie^i)] , (64) 



where 



V_ 
W 



(65) 



We chose to characterize the ring-down time scale as 



tRD = Aq ^ + ^ = 



2W' 



(66) 



Note that in the limit C ~* 1 slowing down occurs and 
tRD — oo. This limit corresponds to the case of operating 
the resonator near a jump point, close to the edge of the 
bistability region. On the other hand, the limit C = 
corresponds to the linear case, for which t^iu at resonance 
{ujp = ujq) is given by Eq. 



VIII. DISPLACEMENT SENSITIVITY 



3V3 ^\K\ - V373) 
and the resonator mode amplitude is 

' ^/3(|X|-^/373) 



c{t)= / dt'Git-t')r{t') 



where 



r (t) = + W*F it) - VF^ (t) 

at 



The propagator is given by 
G{t) = u (t) 



Ai — Ao 



(57) 



(58) 



D. Ring-Down Time 

The solution of the equation of motion (|50p can be 
expressed as [l!] 



(59) 



(61) 



where u{t) is the unit step function, and the Lyapunov 
exponents Aq and Ai are the eigenvalues of the homoge- 
neous equation, which satisfy 



Ao + Ai = 2W' 



(62) 



Consider a measurement in which X^^^ (t) is moni- 
tored in the time interval [0, t] , and the average measured 
value is used to estimate the displacement x. Assuming 
the case where r is much longer than the characteristic 
correlation time of AT^^o (^)j one finds that the minimum 
detectable displacement is given by [l| 



5x = 



dx 



(67) 



Moreover, using Eq. (j43[) . one finds the minimum de- 
tectable displacement per square root of bandwidth, 

1/2 

which is defined as Px — {t/2tt) Sx, is given by 



Prr. 



dlogM 



da; 



R-'P'J' (0) , 



(68) 



where R is the responsivity with respect to a change in 
Wo, namely R — |9Ao/i9a;o|. To evaluate Px we calcu- 
late below the responsivity R and the spectral density 
(60) P^^^ (0) using results, which were obtained in Ref. [2^ . 



A. Responsivity 

In steady state, when input noise is disregarded, the 
amplitude of the fundamental mode of the stripline res- 
onator is expressed as Be'^'^''*, where the complex num- 
ber B is found by solving Eq. (|49p . and the average 
output signal is expressed as houte"^^^* [bout is in gen- 
eral complex). The amplitude of the output signal in the 
feedline bout is related to the input signal bin and the 
mode amplitude B by an input-output relation (23j 



\o\i = \W\^ -\V[' 



(63) 



'B . 



(69) 
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Consider a small change in the resonance frequency 
Suq. The resultant change in the steady state mode am- 
plitude SB can be calculated using Eq. (|49|) 



i (dujo) B + W5B + V (SB)* = 



(70) 



The solution of this equation together with Eqs. and 
([69)1 allows calculating the responsivity 



|sin ((/.t + (/'c) + C sin (0t-(/.c) I 



where 



0c 



20LO - 0w + - 201 



(71) 

(72) 
(73) 



gi</.s ^ 5/ g*0w = 14^/ and e*"^^ = \V\. 



B. Spectral Density 

The zero frequency spectral density Px (0) was calcu- 
lated in Ref. [23| 



PxiO) 



271 (1 - Ce^*"**) 



72 

71 



271 (1 - (e^^-^'' 



coth 



hujQ 

2kBTi 



\w\{i-e 



coth 



hujQ 

2fcsT2 



273 
71 



271 (1 - Ce^'"^') 



IW'I(I-C^) 



coth 



2kBT3 



(74) 



C. The Linear Case 



In this case K — — 0, W = i (wq — Wp) -I- 7 and 
F = 0, and consequently C = 0. Thus, the responsivity 
is given by 



R 



272^1^1 |sin(0t + 0c) 
\W\ 



(75) 



Moreover, for the case where Ti = T2, Eq. (|7^ becomes 



Px (0) = coth 



2fcsri 



(76) 



The largest responsivity is obtained at resonance, namely 
when ujp = luq, and when the homodyne detec- 
tor measures the phase of oscillations, namely, when 



|sin {(pt + 0c) I = 1- For this case Px obtains its smallest 
value, which is denoted as Pxo, and is given by 



PxO 











ax 





7 



871 



■ coth 



^0 



1/2 



(77) 

Note that this result coincides with the SLL value P^^^ 
given by Eq. ([T]) provided that the quality factor in Eq. 
([T]) is taken to be given by Q = 2^071/7^. That is, 
PxO ~ Px^^ in the limit of strongly overcoupled res- 
onator, namely when the damping of the resonator is 
dominated by the coupling to the feedline (71 ~ 7). 

The response of the stripline resonator is approxi- 
mately linear only when |-B| is much smaller than the crit- 
ical value corresponding to the onset of nonlinear bista- 
bility. Using this critical value \B\^, which is given by 
Eq. ((511), and using Eqs. (g^l) and (gi]), one finds that 
the smallest possible value of Pxo in this regime is roughly 
given by 



0.14 



d log A/ 

da; 



i_ 



V373 
li^l 



coth 



2kBTi 
(78) 



D. The General Case 



The minimum detectable displacement per square root 
of bandwidth in the general case can be written as 



Px=Pl 



SLL, 



(79) 



where 



ff(0t) 



M 


■(i-C=)Vx(o)" 


1/2 
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|sin(0t -1- 


0c) + Csin((?! 


t - 0c)| 



(80) 



and Px (0) is given by Eq. ([74|) . The function g{4>t) is 
periodic with a period tt. Its minimum value is denoted 
as f^min- Beating the SLL given by Eq. ([T]) is achieved 
when (7niin is made smaller than unity. 



IX. VALIDITY OF APPROXIMATIONS 

In this section we examine the conditions, which are 
required to justify the approximations made, and deter- 
mine the range of validity of our results. Clearly, our 
analysis breaks down if the driving amplitude bin is made 
sufhciently large. In this case both the adiabatic approxi- 
mation and the assumption that back-reaction effects are 
negligibly small will become invalid. The range of validity 
of the adiabatic approximation is examined below by esti- 
mating the rate of Zener transitions. Moreover, we study 
below the conditions under which back-reaction effects 
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acting back on the mechanical resonator play an impor- 
tant role. Using these results we derive conditions for the 
validity of the above mentioned approximations. These 
conditions are then examined for the case where the sys- 
tem is driven to the onset of nonlinear bistability. This 
analysis allows us to determine whether the device can be 
operated in the regime of nonlinear bistability, where the 
effects of bifurcation amplification 24, 25, 2(3^7 ^8. i29ll 
and noise squeezing can be exploited [23, [3fl, [3l|, |33 |. 
without, however, violating the adiabatic approximation 
and without inducing strong back- reaction effects. 



1 Aoj 

2 Ulm 



2 / ^ \ 2 



1 + ^ 



7 Aix 



2 / \ 2 

J_ ] + f 1 > A, 



(i) 



1 - 



(86a) 



2 / \ 2 



(86b) 



Adiabatic Condition 



and 



As before, consider the case where the externally ap- 
plied flux is given by <i>e = *i'o/2, and the stripline res- 
onator is driven close to the onset of nonlinear bistabil- 
ity, where the number of photons approaches the critical 
value given by Eq. (j58[) . and assume for simplicity the 
case where 73 <C \K\. Using Refs. one finds that 

adiabaticity holds, namely Zener transitions between adi- 
abatic states are unlikely, provided that 



> 1 



where 



27rMo 2haje\B\ 



(81) 



(82) 



(87) 



We now wish to examine whether backreaction effects 
are important when the stripline resonator is driven to 
the onset of nonlinear bistability. For simplicity we ne- 
glect nonlinear damping, that is we take 73 = 0. Using 
Eqs. dUl), (gH) and one finds that the value of Rq 
[Eq. (|85p ] at the onset of nonlinear bistability, which is 
denoted as i?oc, is given by 



Oc 



64A(^)' 



mujt, 



(88) 



It is straightforward to show that |xf| < uj„i/j and |xd| < 
^m/7 for any value of Alu. Thus, back-reaction can be 
considered as negligible when 



B. Back Reaction 



77" 



-Rqc < 1 



(89) 



Back reaction of the driven stripline resonator results 
in a force noise acting on the mechanical resonator and 
a renormalization of the mechanical resonance frequency 
uj,n and the damping rate 7m [3^ . [35| . The renormalized 
values, denoted as uj'^ and 7^, are expressed using the 
renormalization factors R{ and i?d 

^cff ^ 

i?f = ^ , (83a) 

7m 

i?d = ~ , (83b) 

7m 

which where calculated in Ref. [3] 

Ri = XiQmRo , (84a) 
Rd = XdQmRo , (84b) 

where Qm = UJrn/jm, 

8hni\B\U^f 

Ro = 1 — , (85) 



X. EXAMPLE 

We examine below an example of a device having the 
following parameters: Zt = 5017, L}j/LtI — 1, a;e/27r = 

5 GHz, ujeh = lO"*, 72 = O.I71, and 73 = O.liT/VS 
for the stripline resonator, m = 10~"'^^kg, Wm/ZTr = 
0.01 GHz, and LOmllm = 10'^ for the nanomechanical 
resonator, A = 9.1 x lO'i'^H, Cj = 0.86 x lO^^^'F, 
and Ic = 0.44 ^A for the RF SQUID, couphng terms 
Afo/Lfc = 0.005, Lfc/A = 0.25, |dlogA//da;p^ = lOOnm, 
and temperature T — 0.02 K. These parameters for the 
stripline resonator [3^ , for the nanomechanical resonator 
[33|, and for the RF SQUID [H IH, are within reach 
with present day technology. 

Using FastHenry simulation program 
(www.fastfieldsolvers.com) we find that the chosen 
value of A corresponds, for example, to a loop made of 
Al (penetration depth 50 nm) , having a square shape 
with edge length of 190 /xm, line width of 1 /im and film 
thickness of 50 nm. The values of Cj and Ic correspond 
to a junction having a plasma frequency of about 
25 GHz. 
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Using these values one finds (3l = 1-2, f3c = 0.1, and 
Eq/H = 560 GHz. The values of /3i, Pc and um are em- 
ployed for calculating numerically the eigenstates of Eq. 

14]. From these results one finds the parameters 
in the two-level approximation of Hamiltonian Tii [Eq. 
([37]) ] 77 = 1.2Eo and A ^ 0.073£'o, and the nonhnear 
terms ^2/211 = 1.3 MHz and K/2Tr = 0.19kHz. 

The validity of the adiabatic approximation is con- 
firmed by evaluating the factor in inequality (|8ip 



rihVc 



12 



(90) 



whereas, to confirm that back-reaction effects can be con- 
sidered as negligible we evaluate the factor in inequality 



77" 



-Rqc = 1.4 X 10"'' 



(91) 



Using Eq. ([78| we calculate the smallest possible value 
of F.o 



Pm.c = 6.9 X 10 



-14 



Hz 



(92) 



A significant sensitivity enhancement, however, can be 
achieved by exploiting nonlinearity. Figure ([3]) shows 
the dependence of the mode amplitude \B\, the factor 
^, the reflection coefficient \bout/bin\^ , the enhancement 
factor (7min, and the ring-down time ^rd on the pump 
frequency. Panel (a) represents the nearly linear case for 
which bin = 0.001 (6m)c> panel (b) the critical case for 
which bin = (^Mi)c' and panel (c) the over-critical case 
for which 6i„ = 1.5 (6j„)^. As can be clearly seen from 
Fig. ([3]) a significant sensitivity enhancement gmin <SC 1 
can be achieved when driving the resonator close to a 
jump point. However, in the same region, slowing down 
occurs, resulting in a long ring-down time iriD 3> 2/7. 



Ql 1 0' ' 0' ' 
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2 
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FIG. 3: (Color online) Dependence of the mode amplitude 
\B\, the factor C^, the reflection coefficient \bout/bi„f , the en- 
hancement factor gmin, and the ring-down time tao on the 
pump frequency. Panel (a) represents the nearly linear case 
for which bin ~ 0.001 (bin)^, panel (b) the critical case for 
which bin ~ (bin)^, and panel (c) the over-critical case for 
which bin = 1.5 {bi„)^. 



XI. DISCUSSION AND CONCLUSIONS 

In the present paper we consider a resonant detection 
configuration, in which the response of the resonator is 
measured by monitoring an output signal in the feedline, 
which is reflected off the resonator. We show that op- 
erating in the regime of nonlinear response may allow 
a significant enhancement in the sensitivity by driving 
the stripline resonator near a jump point at the edge 
of the region of bistability. The factor gmin, which repre- 
sents this enhancement, can be made significantly smaller 
than unity in this limit. However, this result is not gen- 
eral to all resonance detectors. Consider an alternative 
configuration, in which, instead of measuring an output 
signal reflected off the resonator, the response is mea- 
sured by directly homodyning the mode amplitude in 
the resonator. In the latter case, a similar analysis yields 
that the enhancement factor gmin cannot be made smaller 



than 0.5 [40]. Thus, to take full advantage of nonlinearity, 
it is advisable to monitor the response of the resonator 
by measuring a reflected off signal, rather than measuring 
an internal cavity signal. 

As was discussed above, the largest sensitivity en- 
hancement is obtained close to the edge of the region 
of bistability, where 1 — C ^ 1- On the other hand, in 
the very same region, the perturbative approach, which 
we employ to study the response of the stripline res- 
onator, becomes invalid since the fluctuation around 
steady state, which is assumed to be small, is strongly 
enhanced due to bifurcation amplification of input noise 
[2^ . The integrated spectral density of the fluctuation, 
which was calculated in Eq. (80) of Ref. |4^] (for the 
case Ti = T2 = T3), can be employed to determine the 
range of validity of the perturbative approach 



1 



1-C 



■ coth 



^0 
2kBTi 



« \B\ 



(93) 
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By applying this condition to the case of operating at the 
onset of bistabihty point [see panel (b) of Fig. ^] , one 
finds that the smallest value of ^min in the region where 
inequality holds is ~ lO^'^ for the set of parameters 
chosen in the above considered example . Thus a signif- 
icant sensitivity enhancement is achieved even when the 
region close to the jump, where the perturbative anal- 
ysis breaks down, is excluded. On the other hand, the 
perturbative analysis cannot answer the question what 
is the smallest possible value of gmin- Further work is 
required in order to answer this question and to properly 
describe the system very close to the edge of the region 
of bistability, where nonlinear terms of higher orders be- 
come important. 

In the present paper only the case where back-reaction 
effects do not play an important role is considered. 
Thus our results are inapplicable when the displacement 
sensitivity approaches the value corresponding to zero 
point motion of the mechanical resonator y^h/irujJmlm 



(1.6 X 10~^^ m /VTlz for the above considered example), 
and con seq uently quantum back-reaction becomes im- 
portant [4l|. Moreover, the results are valid only when 
inequality ([M)) is satisfied. To increase the range of appli- 
cability of the theory and to study back-reaction effects, 
a more genera l ap proach will be considered in a forth- 
coming paper [42|. 
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